No-scalar hair conjecture in asymptotic de-Sitter spacetime* 
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We discuss the no-hair conjecture in the presence of a cosmological constant. For the first step 
the real scalar field is considered as the matter field and the spacetime is assumed to be static 
spherically symmetric. If the scalar field is massless or has a convex potential such as a mass term, 
it is proved that there is no regular black hole solution. For a general positive potential, we search 
for black hole solutions which support the scalar field with a double well potential, and find them 
by numerical calculations. The existence of such solutions depends on the values of the vacuum 
QQ ' expectation value and the self-coupling constant of the scalar field. When we take the zero horizon 

Q\ , radius limit, the solution becomes a boson star like solution which we found before. However new 

ON ■ solutions are found to be unstable against the linear perturbation. As a result we can conclude that 

the no-scalar hair conjecture holds in the case of scalar fields with a convex or double well potential. 

Or 

ON . 

I. INTRODUCTION 

>. 

' Recent developments of observational techniques have increased interest about black holes and increasing data 
shows evidence of supermassive black holes in the center of galaxies, and solar mass size black holes which form 
binary systems. From the theoretical view point, however, there remain many unsolved fundamental issues related 
to the black holes. One of them is the validity of the black hole no-hair conjecture proposed by Ruffini and Wheeler 
0. The black hole no-hair conjecture states that after the gravitational collapse of the matter field, the resultant 
[ black hole approaches stationary spacetime and that all its multipole moments are then uniquely determined by two 
parameters, M and a, which are physically interpreted as the mass and the angular momentum of the black hole. 
When the source has a net charge Q, then of course the parameter is also requested to uniquely determine its (electric 
and gravitational) multipole moments. The black hole no-hair conjecture is supported by the black hole uniqueness 
£jpj| theorems in electrovacuum theories and by the works of Chase ||, Bekenstein Q, Hartle || and Teitelboim p 
which show that stationary black hole solutions are hairless in a variety of theories coupling classical fields to Einstein 
^ ■ gravity. 

However a non-trivial static, spherically symmetric solution, called colored black hole, was discovered in the Einstein- 
Yang- Mills system . It can be interpreted as the self-gravity of the Yang-Mills field being supported by its repulsive 
force. Although this solution is found to be unstable both in the gravitational sector B and in the sphaleron sector [p| , 
non-Abelian hair is generic and many other non-Abelian black holes were discovered after the colored black hole. Note 
that some of them are stable and that such solutions are the counterexamples of the black hole no-hair conjecture. The 
essential difference between the field equations of Einstein- Yang-Mills systems and those of the Einstein-Klein- Goldon 
system is in the form of the potential term. By using these properties, Bekenstein J|,|lO| and Sudarsky Jll]] provided 
simple proofs for the no-scalar hair theorem in spherically symmetric asymptotically flat spacetime, in the case where 
the matter consists of a single scalar field with a convex potential and in the extended case where the matter consists 
of multiple scalar fields with an arbitrary positive potential. Heusler also proved the no-scalar hair theorem by using 
a scaling technique [^2[ . 

Most of the proofs of the no-hair theorems impose, however, asymptotic flatness. Hence the following natural 
question arises: Can we extend no-hair theorems to spacetimes with different asymptotic structure? This is the main 
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issue of this paper. Here we consider the system including the cosmological constant and only one real scalar field 
as a matter field for the first step. The importance of the cosmological constant comes not only from the theoretical 
aspect but also from the observational results of our universe. For example some astrophysicists have pointed out that 
the small cosmological constant may explain the observed number count of galaxies |13| , and the recent observation 
of the type la supernova at large redshift also supports the universe with a cosmological constant [O]. Furthermore, 
in the early universe, we usually expect a vacuum energy, which is equivalent to the cosmological constant. As for the 
black hole solutions with a cosmological constant, the family of the Kerr-Newman-de Sitter solutions are known as the 
exact solutions. Besides these the cosmic colored black hole solution is derived |HJ. These solutions are an extension 
of the Kerr-Newman solutions and the colored black hole to the non-zero cosmological constant case. Interestingly, 
although there are black hole solutions in the Einstein- Maxwell-dilaton system |L6[ , no black hole solution exists if 
we take the cosmological constant into account [|l7j . Hence the cosmological constant can strongly affect the existence 
of black hole solutions. 

If we find the physical solutions with scalar hair, the no-scalar hair conjecture can not hold in asymptotic de Sitter 
spacetime and such solutions may have an influence on the cosmology. If there are no such solutions, we can say that 
the no-scalar hair conjecture holds in our system. In other words, the weak cosmic no-hair conjecture holds in the 
sense that all initially expanding universes with positive cosmological constant approach the "de Sitter spacetime" 
asymptotically, except for locally distributed black holes. Hence it must be stressed that the black hole no-hair 
conjecture in asymptotic de Sitter spacetime is strongly related to the weak cosmic no-hair conjecture. 

This paper is organized as follows. In section II we introduce the model and the basic equations and try to extend 
the no-hair theorem in a spherically symmetric system. In section III we derive black hole solutions numerically in the 
model that the real scalar field has a double well potential and discuss their properties. We investigate the stability 
of these solutions in section IV. We give our conclusions and remarks in the final section. 



II. MODEL AND BASIC EQUATIONS 



We will consider the model given by the action 

S = I d^x\]—g 



^(*-2A)-±(V*)'-V(0 



(1) 



where 4> is the real scalar field and V((f>) is its potential. We shall assume a spherically symmetric spacetime and 
adopt the Schwarzschild type metric, 



ds 2 



1 



2Gr 
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dt 2 
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dr 2 
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The mass function m and the lapse function 5 depend on both the time coordinate t and the radial coordinate r. 
The mass function is the quasilocal mass defined by Ref. Jl9(| , which is the gravitational energy subtracted by the 
energy due to the cosmological constant A. In other words, it is considered as the energy of the matter field. In 
spherically symmetric spacetime, the energy m is nondecreasing in the outgoing null or spacelike direction in the 
region 1 — 2Gm/r — 3Ar 2 > if the matter fields satisfy the dominant energy condition. 
Varying the action (Q) and substituting Ansatz (||), we derive the field equations 



Here, we have used the dimensionless variables, t = y/Xt, f = \/~Rr, m = y/AGm, <fi = y/Gcj) and V = GV/A. A dot 
and a prime in the field equations denote derivatives with respect to t and f , respectively. In this section and the 
next section we consider only the static solutions, hence we drop the time derivative term of the field equations for a 
while. 

For the boundary conditions of the metric functions on the black hole event horizon (BEH) and on the cosmological 
event horizon (CEH), we impose the following three Ansdtze. 
(i) The existence of a regular BEH tb] i.e., 



2Gm{r B )=r B - 
5{tb) < oo. 
(ii) The existence of a regular CEH rc; i.e., 

2Gm(r c ) = 2GM = r c (l 
S(rc) < oo. 



A 



(7) 
(8) 

(9) 
(10) 



Here we assume 5(rc) — 0. Note that if we are interested in a different boundary condition such as S — » 5* ^ 0, we 
can always have such a boundary condition without further calculation by rescaling the time coordinate. That is, 
introducing 5 = 8 — 5* , and rescaling the time coordinate as t = e~ s t, we recover our boundary condition, 
(iii) The non-existence of singularity between BEH and CEH; i.e., for Tb < r < rc> 



A .- 

2Gm(r) < r I 1 - —r" 



(11) 



As for the scalar field, we impose the finiteness of itself and its derivatives. These conditions guarantee that the 
curvature invariant I = RfivpaR^ 1 " 7 is finite and that no naked singularity appears. 

Now we investigate whether or not the no-scalar hair theorems obtained by Bekenstein jl 10 and Sudarsky (ll) can 
be extended to our system. First we consider the case of the scalar field with convex potential, e.g. V((f>) — m^cf) 2 . In 
the static spacetime the equation of the scalar field is 



(f 2 e- 5 f4>')' 



r 2 e~ 5 



dV(4>) 



(12) 



where / = 1 — 2m/f — f 2 /3. Multiplying by <j> and integrating from BEH to CEH, we obtain 



dr 



rc 

tb 
rc 

rB 



dff 2 e 



drr 2 e 



-~^{f 2 e- s f^)' 

'~JV{4>) 



~dv(4>) 



0'- 



r 2 e~ 5 



rc 



(13) 



Since V is a convex potential, the integrand is positive semi-definite. Hence the possible solution is <j) = 0. As a 
result the same argument as Bekenstein's holds and the no-scalar hair theorem is complete in the convex potential 
case H(|. When the scalar field is massless, i.e., V = 0, Eq. ( |l3| ) also demands = 0, which is the no-hair theorem 
for a massless scalar field. 

Next we proceed to the theorem proved by Sudarsky. First we rewrite Eq. (^) in the static spacetime as 



Multiplying <f)' and integrating once, we obtain 



e-^ = 0. 



(14) 



E' = -ae 



-<5 1/2 



(15) 
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where 



E = e- s (/0' 2 - V) , (16) 
2 / 3m 1 ,\ 

In the A = case, i.e., asymptotically flat case, E — —V < on BEH, and a is always positive between BEH and 
CEH because of the regularity. (Note that we have normalized variables by A, a — 2(1 — 3m/2r)/r in the A = 
case.) Hence E decreases monotonically and then becomes negative everywhere. However E must approach zero in 
the asymptotic region JTlJ] . This means that there is no regular solution in this system. On the other hand, in the 
A 7^ case, E = —V < on both horizons and a is no longer positive definite. Hence the inconsistent behavior of E 
can not be derived from Eq. (|l5|). Furthermore asymptotic behavior does not reject the scalar hair in the arbitrarily 
positive potential case 

This suggests that the method adopted here may not be appropriate to prove the black hole no-hair or that the 
black hole solutions with non-trivial scalar hairs may exist. We will search for such solutions in the next section by 
using numerical calculations. 



III. BLACK HOLE SOLUTIONS WITH SCALAR HAIR 



In this section we investigate the system including a double well potential V(4>) — A(cf) 2 — v 2 ) 2 /4 as an example 
of a positive potential. In this case there are two trivial black hole solutions. One is (a) <fi = v, m = M, S = 0. 
The scalar field takes its vacuum expectation value and does not contribute to the spacetime at all. This is the 
usual Schwarzschild-de Sitter solution. The other is (b) cf> = 0, m = M + TrXv 2 f 3 /3, (5 = 0, where A = A/GA and 
v = \/~Gv. The scalar field sits on the top of the potential barrier everywhere and its contribution to the spacetime 
can be interpreted as the "effective cosmological constant" defined by A e /y = 2ir\v 4 . We call this solution the excited 
Schwarzschild-de Sitter solution. 

The equation of the scalar field (|]) is rewritten as 

4>' = \4> U 2 - s 2 ) . (18) 

On the event horizon Eq. ( fl8| ) becomes 

f'4>' = \4> (4? - i 2 ) . (19) 

Since /' is positive on the BEH, </f>'(r\B) > if — v < </>(rs) < 0, v < 0(rs), and </>'(rs) < if <j)(rB) < —v, 
< <f)(rB) < v- On the other hand /' is negative on the CEH, 0'(rc) > if <j)(rc) < —v, < (j)(rc) < v and 
4>'{ r c) < if — v < 4>{ r c) < 0, v < 4>{rc)- At the extremum point of 0, Eq. ( [i"8| ) becomes 

f4>" = \~4> (ft - v 2 ) . (20) 

Since / is always positive between BEH and CEH, <f> has no local maximum if — v < <fi(r) < 0, v < (j)(r) and <f) 
has no local minimum if 4>{r) < —v, < cf>(r) < v. Now we estimate the possible behavior of the scalar field from 
the above constraints. We can restrict 4>(rB) > without loss of generality because the scalar field has reflection 
symmetry. If 4>{ t b) > v, </> increases around the BEH. Since there is no maximum point cf> > v, it continues to increases 
monotonically to the CEH. However on the CEH <fi' must be negative. This is a contradiction. When < </>(rs) < v, 
4> decreases around the BEH. Since there is no minimum in this region, 4> must pass over the potential barrier (0 = 0). 
Continuing to decrease beyond <j) — —v, <j) can not satisfy the boundary condition on the CEH because <j) has no 
minimum in <fi < —v. Hence <j) must stop at a certain value in — v < <fi < or have a minimum. In the latter case <p 
passes over the potential barrier again. In this way <j> may oscillate any number of times and takes a value —v<<fi<v 
on the CEH. For each oscillation <p has to go over the potential barrier. 

Now we search for non-trivial static solutions using numerical analysis. We drop the time derivative term of the 
field equations (||), (Q) and (||) and integrate them from the BEH with the boundary conditions (^) and (|8|). Since 
the equation of the scalar field (|^) becomes singular on the event horizons, we expand the equation and variables by 
power series of r — rs to guarantee the regularity on the BEH, and use their analytic solutions for the first step of 



f<t>" 



- + 4irrif>' 2 
f 



f + r 
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integration. For most of the values of <^(fs), the scalar field diverges as the integration approaches the CEH. Hence 
we have to find a suitable value of ^>{tb) m order to satisfy the boundary conditions ([)]) and (10). In this sense 4>(rB) 
is a shooting parameter. 

We found the non-trivial solutions when A, v and ?b satisfy a certain condition which will discussed later. They 
are classified into several families by the node number n of the scalar field. Configurations of the field functions of 
the solutions with A = 700, v = 0.1 and n = 1 for different radii of BEH are shown in Fig. 1. The left and right end 
points of each solid line are BEH and CEH, respectively. A dotted line traces them. A dashed line is the regular 
solution, which can be considered as the solutions in vb —* limit plj . The structures of the new solutions do not 
concentrate around the BEH but spread out to a cosmological scale!~As the BEH becomes large, the CEH becomes 
small and the values of the scalar field on each event horizon approaches its vacuum value <p = ±5. Hence the scalar 
field must vary rapidly in a small region and become steep. However, as the BEH becomes even larger, </>(?\b) and 
4>{fc) approach zero and finally BEH and CEH coincide to be the extremal solution. 

We show the A dependence of the scalar field for ?b — 0.2, v — 0.1 and n = 1 in Fig. 2. For large A, which can be 
considered as the case that the cosmological constant is small, the scalar field remains at its vacuum value <j) & v even 
for quite large f around the BEH. As A decreases, the amplitude of <j) becomes small and finally the solution coincides 
with the excited Schwarzschild-de Sitter solution (b) at a certain critical value X m in — 354.8. These properties are 
qualitatively the same as the regular solution discussed in Ref. |pif . For different values of parameters v, we obtain 
similar results to those mentioned above. 

We show the M — fs diagram of new solutions in Fig. 3, where M is the quasi-local mass on the CEH. We fix 
v = 0.1 and n — 1. We also show (a) the Schwarzschild-de Sitter branch and (b) the exited Schwarzschild-de Sitter 
branch by dashed and dotted lines, respectively, for comparison. We find that each solution branch turns at M = 1/3. 
From the boundary condition ([)]) 

Since the BEH has a one to one correspondence to the CEH around the turning point, 

or B or B or c 2ar B 

on the turning point. From Fig. 3, drc/dr B does not vanish, hence fc — \fhrc = 1 at the turning point. This means 
M = 1/3. 

The properties of the new solution branch are quite different depending on A. For small A there is only one solution 
for each BEH radius. As the BEH radius increases, the branch approaches the excited Schwarzschild-de Sitter branch 
and finally both branches coincide at the maximum BEH radius. At this point the BEH and CEH coincide and the 
solution becomes extremal. There is the minimum BEH radius, below which only trivial solutions, i.e., Schwarzschild- 
de Sitter solutions (a) and (b) exist. The non-trivial branch is very close to the excited Schwarzschild-de Sitter 
branch. As we can see in Fig. 2, the scalar field remains near the top of the potential barrier (<fi ~ 0) for small A. The 
main contribution to M must be this potential energy of the scalar field. Hence the solution is similar to the excited 
Schwarzschild-de Sitter solution. 

On the other hand, for large A, there are two solutions, which have the same BEH radius. This only occurs for a 
narrow range of radii The maximum BEH radius is not at the extremal point but a turning point plotted by dots in 
Fig. 3. On the small BEH radius , we can take the rs — * limit even for the non-trivial branch, and the solution 
becomes the regular solution pl|. From Fig. 2 the scalar field takes an almost vacuum value near both horizons and 
passes rapidly on the top of the potential barrier. As a result the non-trivial branch is quite different from the excited 
Schwarzschild-de Sitter solution. 

In Fig. 4 we show the critical parameters in the A — v plane for rg = 0.6 and n = 1. There is a solution only 
in the region between X m im and X max . The existence of the critical line X m in plotted by the dashed line is due to 
reasons similar to those in the regular solution case. In the regular solution case the critical lines are derived by 
investigating the perturbation from the de Sitter solution and are expressed as (R CO s/ Xcomp) 2 = n(2n + 3), where 



Rcos — Y 3/(1 + 27rAw 4 ) is the typical scale of the CEH, Xcomp = y2/Xv 2 is the typical size of the structure and 
n is the node number of the scalar field [£lj . If A becomes small, the size of the structure becomes large compared 
to the size of the CEH. Below the critical parameter X cr the structure can not be packed into the radius R coa and 
the non-trivial solution disappears. Let us consider the effect of the existence of the BEH to the critical parameter. 
Rescaling f and rh by the BEH radius fs as f = f/rs and m = fh/friB, and taking an approximation 3> </>' ~ 0, 
i.e., near the critical parameter A m j„, the field equations become 
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^ m ~ \~2 -2/12 ~2\2 / O o\ 

— w Ar B r (0 - u ) , (23) 

Af|0(^ 2 - « 2 ) « 0, (25) 

where Af 2 ^ can be considered as the effective self-couphng constant and the critical solutions are controlled by it. 
Consequently \ m in becomes small when we consider the large BEH radius. 

On the other hand X max is a new restriction which can not be seen in the regular solution case. From Fig. 3 the 
maximum BEH radius for each A branch becomes small as A increases. Hence if we fix the BEH radius, the non-trivial 
solution disappears at a certain critical value of A. This is a X max curve. The types of the critical solutions are 
different depending on X rnax or equivalently v max = v(X = X max )- When v ma x^p.l4: for Tb — 0.6 the critical solution 
is the extremal solution. Otherwise the critical solution is not the extremal one, but is shown by the points plotted 
on the solution for A = 2000, 5000 in Fig. 3. 



IV. STABILITY ANALYSIS 



In the previous section we found black hole solutions with scalar hair in de Sitter spacetime although there is no 
counterpart in asymptotically flat spacetime. This means that the no-hair conjecture may not hold in asymptotic de 
Sitter spacetime. In this section we investigate the stability of new solutions by using a linear perturbation method 
in order to check whether the scalar hair is really physical or not. Here we focus only on the radial modes. 

First we expand the field functions around the static solution <pQ, mo and Sq as follows: 



r 

m(t, f) = rho(t) + mi(t, r)e, 
5(t,f) = S (t) + 5 1 (t,r)e. 

Here e is an infinitesimal parameter. Substituting them into the field functions 
and higher order terms of e, we find 

mi = 47rf 2 /o^ ^i! 



(26) 

(27) 
(28) 

(||) and dropping the second 
(29) 



e*7o~Vi + 



-So 



Ml 



r (e- Ao /o) + Snre^X^l - « 2 )0 O O + Xe^^l - v 2 ) 



\ (re-toft)' 



fhi = 0, 



(30) 



where /o = 1 — 2rho/f — 1/3. Next we set 4>\ = £(f)e zcrt and fh\ = i](r)e i ' 7t . If a is real, <j> oscillates around the static 
solution and then the solution is stable. On the other hand, if the imaginary part of a is negative, the perturbation 
(j>i and tiii diverges exponentially with time and then the solution is unstable. By Eq. (|29| ) the relation between £ 
and rj is rj = 47rr/b0o£- Then the perturbation equation of the scalar field becomes 



d 2 t 



df* 2 

where we employ the tortoise coordinate f* defined by 

dr* 
df 



/o' 



(31) 



(32) 



and the potential function is 
U(r) = 



-So 



fa 



-So 



fo) +87 



X^l-v^o^ + Xe-toiWl-v 2 ) 



(33) 
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Fig. 5(a) shows the potential functions U(r) of the solution with A = 300, v = 0.1 and n = 1. 

Since d 2 ^/df* 2 = U(r) = on both horizons, £ must approach zero as f* — > ±oo for the negative mode by the 



regularity of Eq. (31). Under this boundary condition we have searched for the negative eigenmodes and found them 
as shown in Fig. 5(b). These modes are bound states (m = 0) and there is no excited mode (m = 1) for n = 1. For 
the solution with n = 2, however, there are two negative modes. We expect that the node number of the scalar field 
of the static solution n exactly corresponds to the number of their negative modes m. Fig. 5(c) gives the eigenvalues. 
We also plot the eigenvalues of the excited Schwarzschild-de Sitter case with dashed lines for comparison. At the 
point P, where the new solution branch coincides with the Schwarzschild-de Sitter branch, another unstable mode, 
which is an excited mode, appears in the Schwarzschild-de Sitter branch. This result is consistent with the analysis 
using catastrophe theory P4P3| , and this stability change is classified into the swallow tail catastrophe. Although 
there is the tb — > limit for large A as we mentioned, the eigenvalue does not vary continuously to that of the regular 
solution because of the difference of the boundary condition around the origin. Varying the parameters we found 
negative eigenmodes for every non-trivial solution. As a result, all of the new solutions are unstable, even against the 
radial perturbations and the scalar hair falls out easily. Thus, although we found black hole solutions with scalar hair 
in the presence of cosmological constant, such hair is not physical and the no-hair conjecture seems to hold even in 
the de Sitter spacetime. 

V. CONCLUSION 

We examined the no-hair conjecture in the presence of a cosmological constant. For the first step, the real scalar 
field was considered as the matter field and the spacetime was assumed to be static spherically symmetric. When 
the scalar field is massless or has a convex potential like mass term, it was proved that there is no regular black hole 
solution. However we can not find any proof excluding scalar hair in the general positive potential case. Therefore we 
searched for black hole solutions which have a scalar field with a double well potential, and found them by numerical 
calculations. Their field configurations spread out to the cosmological scale and are classified by the node number of 
the scalar field. For large BEH a solution branch ends up with the extreme solution in the M-fs diagram, while for 
small BEH the behavior is different, depending on the self- coupling constant A. If A is large, we can take the limit 
Tb ~ > and the corresponding solution becomes regular a solution without the BEH. If A is small, the new solution 
branch hits the excited Schwarzschild-de Sitter branch at a non-zero BEH radius. The new solutions have critical 
parameters \ m in and X ma x between which there are non-trivial solutions. X m in is determined by the ratio of the size 
of the structure to the size of the CEH, while X max comes from the extremes of the solutions. 

In order to specify whether the scalar hair we found is physical or not, we investigated the stability of new solutions 
by using a linear perturbation method. As a result all of the new solutions have negative eigenmodes and were found 
to be unstable. Thus the scalar hair is not real but a wig which falls off easily. Although we have considered only one 
real scalar field with a double well potential, we expect that the general no-scalar hair conjecture holds even if the 
cosmological constant exists. 

Since these new solutions are unstable, the scalar field will be swallowed by the black hole and/or will escape to 
infinity over CEH. Then the solution becomes the stable Schwarzschild-de Sitter solution (a) with the same or smaller 
mass than the initial one. It is interesting to consider the development of the solution with maximum mass M = 1/3. 
Although a part of the scalar field escapes to infinity in the general case, it would be possible to set up the initial data 
where all the energy of the scalar field collapses into the black hole. If the third law of black hole thermodynamics 
is valid even in the present case, it should take an infinite amount of time for the black hole to swallow the scalar 
field and for the BEH to become degenerate with the CEH. It has been shown that similar phenomena occur in the 
evolution of the Kastro-Traschen solution (^4|, which can be interpreted as black holes with Q = M are balanced to 
each other in CEH. For small black holes, they can coalesce to form a larger black hole. However when the size of the 
black holes is larger than a certain critical radius, coalescence does not occur even if we set the initial data to give 
the black holes a large initial velocity toward each other J2l| . Our case would correspond to the critical case that the 
resultant BEH radius coincides with CEH. 

As for the rotating case, the no-scalar hair conjecture can not been proved by using the techniques adopted 



in Ref. j|,[10-12| even in the asymptotically flat case because they strongly depend on spacetime symmetry, i.e., 
staticity and spherical symmetry. On the other hand there is the beautiful result of the uniqueness theorem in the 
asymptotically flat case [|]-[§. ^ seems worth investigating its counterpart in the asymptotic de Sitter case. At a first 
glance, however, we will soon find that the cosmological constant prevents us from constructing Ernst type equations 
and that different approaches are needed. We leave them open questions. 
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Figure Captions 



FIG. 1: The configurations of the scalar field 4> = VG<p with one node (n — 1). Wc set v = \/Gv = 

0.1, A = A/GA = 700 and show the solutions for f B = \fKr B = 0.2, 0.4, 0.6, 0.8. Wc also 
plot the regular solution without BEH by the dashed line for comparison. The left and right 
end points of each lines arc BEH and CEH, respectively, and a dotted line traces them. As 
the BEH radius becomes large, the BEH coincides with the CEH and the solution becomes 
extremal. 

FIG. 2: The configurations of the scalar field (f> — \[G§ with one node (n = 1). Wc set v = \/Gv = 

0.1, r B = \/Ar s = 0.2 and show the solutions for A = A/GA = 400, 700, 2000, 5000. As A 
becomes small, (f> approaches zero and the solution coincides with the excited Schwarzschild- 
de Sitter solution (b) at A = X rn in- 

FIG. 3: The mass-BEH radius diagram of the new solutions. We set v = \^Gv = 0.1 and show the 

solutions for A = A/GA = 300, 400, 700, 2000, 5000. The dashed and dotted lines are the 
Schwarzschild-de Sitter solutions (a) and the excited Schwarzschild-de Sitter solutions (b), 
respectively. For small A there is a lower limit of the BEH, while wc can take the BEH zero 
limit for the solution with large A. The new solution branch and the exited Schwarzschild-de 
Sitter branch approach each other for large BEH and finally coincide. At this point the BEH 
and CEH degenerate and the solution becomes extremal. 

FIG. 4: The critical parameters of the non-trivial solutions with tb = VAYb = 0.6 and n = 1. 

Between X ma x and \ m in curves there exist non-trivial solutions. We also plot the critical 
parameters of the regular solution by the dashed line for comparison. On the critical lines 
A m i m the non- trivial solutions coincide with the excited Schwarzschild-de Sitter solution 
(b) . On the other hand the non-trivial solutions become extremal on the critical lines \ max 
with w>0.14, otherwise they correspond to the solution dotted in Fig. 3. 

FIG. 5: (a) The configurations of the potential function U of the linear perturbation equation for 

v = y/Gv = 0.1, A = A/GA = 300, n = 1 and r B = VAr B = 0.4, 0.6, 0.8. (b) The 
configurations of eigenmodes of the perturbation equation with the same parameters as (a). 
We find the only bound state (m = 1). The eigenvalue for each mode is <r 2 = —0.4165, 
—0.1522, —0.0298, respectively, (c) The eigenvalue of the linear perturbation equation for 
v = \/Gv = 0.1, A = A/GA = 300 and n = 1. Dotted lines are those of the excited 
Schwarzschild-de Sitter solution. At the point P where the non-trivial solution disappears, 
the excited mode of the Schwarzschild-de Sitter solution (m = 1) appears. This is consistent 
with the analysis made using catastrophe theory. 
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Fig. 5(b) 
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